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METHOD AND APPARATU S FOR BACKS CATTER SPECTROSCOPY 



C oooij present invention relates to a method of determining of a physical feature of a 
5 medium, comprising: 

^'producing radiation w ith a light source;^ 

^placing a probe on a sample of rseM. medium, the probe comprising a first optical fiber 

f?y€ 

having a first diameter, and at least a second optical fiber having a secon d diameter; ^ 
^sending light coming from title light source, through the first optical fiber; 
10 ^^collecting first backscattered radiation through the first optical fiber and second 
backscattered radiation through the second optical fiber; 
^producing a first signal based on the first backscattered radiation, and a second signal 
based on the second backscattered radiation; q o c/ ^^f*^^ ^ f ^3 > 
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Such a method is known freaa^metkik ct al^l^ 1-^b» a special device j^used to ^7fi* S^/trfn-J 
determine particle sizes in superficial layers. The device is suitable for measuring ^tjffo\ n ^**** 

particle sizes iij for example,an aqueous suspension with polystyrene spheres, but is not > f)^ f^Jl Qjyhc {s 



fitted to accurately measure particle sizes in living tissue. So, determining whether 
living tissue is normal or precancerous, by way of measuring particle sizes in living 
tissue is not very promising. 
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vivo using a spectroscopic arrangement with ten optical fibers. One of the fibers is 
to irradiate a sample, and nine other fibers collect the reflected light By using a 
multitude of fibers to collect the reflected light, it is possible to calculate scattering and 
absorption coefficients of the sample. However, the method is not suitable for locally 
measuring the optical properties of the tissue. In particular, only mean values of the 
absorption coefficient of a relatively large part of the sample can be determined. 

It is an object of the present invention to locally measure a physical feature, such 
as a concentration, of a substance in a medium. 

The object is achieved by a method as described above, characterized by 
y calculating the physical feature by curve fitting sam measured differential backscatter 
signal to a backscatter function, in which the backscatter function is a function of an 
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average path-length travelled by detected scattered photons, wherein the average path- 
length is independent from an absorption coefficient of the medium, and from a 
scattering coefficient of the medium. Contrary to methods using diffusely scattered 
photons sucn^JDoomBos etaQg]7i5th e method according to the invention, the local 
5 absorption coefficient of the sample is measured in an absolute way, independent of the 
magnitude of the local scattering and absorption coefficients. This facilitates the 
measurement of absolute concentrations of absorbing molecules in a sample without 
requiring prior knowledge of the magnitude of the scattering and absorption 
coefficients of the medium. 
10 £fc<*nj In gu^^Bodiment, the average path-length is proportional to the first fiber 

diameter. This has as additional advantage that the average path-length and thereby the 
average penetration depth into the sample of the photons that contribute to the 

■ 

differential backscatter signal can be controlled by choosing the fiber diameter. As a 

result, the sampling volume can be controlled by adjusting the fiber diameter. 
1 5 Hence, the fiberoptic probe can be engineered to match the relevant dimensions of the 

medium under investigation. 
£ OQ0 ^3 In a particular embodiment, the physical feature is a concentration of at least one 

substance in the medium. 
C a 00 9\ Tte invention also relates to a device for determining a physical feature of a 

20 m edium, c omprising^ 

/ )^a light source for produ cing radiation ^) 

a probe with at least a first and a second optical fiber, the first optical fiber having a 
first diameter and being arranged to deliver the radiation on a sample of s^acT medium 
and to collect first backscattered radiation from^ sample, the second optical fiber 
25 having a second diameter and being arranged to collect second backscattered radiation, 
wherein the second optical fiber is positioned alongside the first o ptical fiber ^> 
^a spectrometer for producing a first signal based on the first backscattered radiation, 
and for producing a second signal based on the second backscattered radiation; 
^processor arranged to determine a measured differential backscatter signal as a 
30 function of wavelength using the first and second signals, 

^ characterized in that the processor is arrangedlo calculate the physical feature by curve 
fitting the measured differential backscatter signal to a backscatter function, in which 
the backscatter function is a function of an average path-length travelled by detected 
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scattered photons, the average path-length being independent from an absorption 
coefficient of the medium, and from a scattering coefficient of the medium. 
£qO(C>2 Furt 

-flanri to a data ca 
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^otQ In another aspect of the invention, the invention relates to a method of^* v°^c#Sr& 

determining a ph ysical featae jrfa mediu m, comprising^ 

^producing radiation with a light source; 

^placing a probe on a sample otffiemeium, the probe comprising a first optical fiber 
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having a first dia meter, and at least a se cond optical fiber having a second diameter; 



'sending light coming from the light source, through the first optic* 
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Collecting first backscattered radiation through the first optical fiber and second 
backscattered radiation through the second optical fiben£> 

/^producing a first signal based on the first backscattered radiation, and a second signal 

Xietefmms a measured differential backscatter signal as a function of wavelength using 
the first and second signals, 
Characterized by 



^N^calculating the physical feature by curve fitting the measured differential backscatter 
signal to a backscatter function, in which the backscatter function is a function of a 
mean free path of photons. In this method, it is assumed that only singly scattered 
photons contribute to the differential backscatter signal and as a result the backscatter 




nt, the physical feature is a concentration of at least one 
25 substance in the medium. ( '^zP^ ? f^ K ffl*2%jT* , M * , 



The invention also relates to a devic^ 
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The present invention will be described below with reference to exemplary 
embodiments and the accompanying schematic drawings, in which: 
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pO* I 73 Fig. 1 is a schematic diagram of a measuring device according to a preferred 

embodiment; 

£ 0 0 / f J Kg. 2a and 2b show cross-sections of a sample and two fiber tips in the situation 
wherein the mean free path of the photons is much larger than the diameter of the 
5 fibers; 

f o Q( 73 Fig- 3 shows the results of Monte Carlo simulations for a homogeneous medium; 
£ q n tVJ Fig. 41a differential backscatter signal normalized at zero absorption for several 

scattering coefficients; 
t<*°1 0 Fig. 5 shows a differential backscatter signal of a dilute suspension of 0.2 \xm 

10 polystyrene spheres^ 
fo&it') Fig. 6 shows the total differential backscatter signal as a function of the reflection 
coefficient fa(X) in the range of 10-100 mm" 1 ; 
£ o o^J^ Fig. 7 shows the measured and calculated average path length rasa function of the 

average scattering coefficient; 
COorvJ 15 Fig. 8 is a graph of measurement results for three different absorption coefficients fa 

showing the dc-fiber signal /, the c-fiber signal J and the differential backscatter signal 
Rbs as a function of wavelength; 
£ O 0 'V*~7 Fig. 9 shows a typical spectrum of an absorption curve A along with the specific 

absorption coefficient of Evans Blue dye; 
Pen "rf>]2Q Fig. 10 shows a measured A * as a function of the absorption coefficient fa at X = 600 

nm; 

r 0 ^l] Fig. 1 1 shows typical spectra measured in a suspension of 1 .0 jxm polystyrene spheres 

with and without Evans Blue dye; 
J2° *^^J Fig. 12 graphically shows a molar extinction coefficient as a function of wavelength; 
jr>t*j25 Fig. 13 shows in vivo measurements and a fit ofthe differential backscatter signal** 

in a human trachea realized using a fiber diameter of 400 prn^ and 
pg * 3 Fig. 14 shows in vivo measurements and a fit of the differential backscatter signal 

in a human trachea showing very low oxygenation indicative for lung tumor. y/ ^ A „ 




£00 3Z] 30 A schematic diagram of a preferred embodiment according to the invention is 

shown in figure 1 . The setup consists of a set of optical fibers for the delivery and 
collection of light to and from a sample 1 under investigation. Light from a light source 
2, for example a Tungsten Halogen lamp (Avantes HL-2000-FHSA), is led through a 

Coo 3 a ti* fy^e^-f fhpr^^ 

AH^W ?CT//YL2L^/ocfrG>?7 , -fit** on Sr/>w&r 2 7, 1^^ t*h<c£tj 



5/02905 



PCT/NL200 I /000657 



10 



15 



20 




first aim 3 of a bifurcated optical fiber. The bifurcated optical fiber is at a distal end 4 
coupled to a first distal end of a delivery- and-collection fiber 5 (in the following 
referred to as dc-fiber 5) which is small enough to be fit through a working channel of a 
clinical endoscope, not shown. A second distal end of the delivery- and-collection fiber 
5 contacts the sample 1. Alongside the dc-fiber 5, a collection fiber 6 is arranged to 
collect light reflected by the sample 1 . The collection fiber 6 (referred to as c-fiber 6) is 
connected to a slave channel of a dual-channel spectrometer 7, for example an Avantes 
SD2000. Preferably, the dc-fiber 5 is polished at a small angle to reduce specular 
reflections. » 

Light reflected back from the sample 1 into the c-fiber 6 is led directly into the 
slave channel of the dual-channel spectrometer 7. A second arm 8 of the bifurcated 
fiber is connected to a master channel of the dual-channel spectrometer 7. Light 
reflected into the dc-fiber 5 is coupled back into the bifurcated fiber, and reaches the 
dual-channel spectrometer 7 via the second arm 8 of the bifurcated fiber. An output of 
the spectrometer 7 is connected to an input of a processor 9 which is arranged to 
analyze signals from the spectrometer 7. 

If only the dc-fiber 5 wotriAfee used to deliver and collect light to and from the 
sample 1, a large fraction of collected light is due to single backscattering from small 
sample depths, sw^rX'singl^to-multiple scattenngratio depends on the scattering 
coefficient and phase function of the sample 1 and on a diameter of the dc-fiber 5. The 
contribution of multiply scattered light to the signal of the dc-fiber 5 can be 
approximately determined by combining the signal of the dc-fiber 5 with a signal 
coming from an additional fiber, i.e. the c-fiber 6 mentioned above. 

a differential backscatter signal R bs as a function of the wavelength X is 
determined using a formula like 



Rbs Q)- c 
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(I(A)-IJX)) 



J(X) 
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where 1(2) is the signal from the dc-fiber 5 in contact with the sample 1, I n (X) is the 
signal from the dc-fiber 5 submersed in a fluid with an appropriate refractive index (for 
tissue: water would be appropriate), 1^ (X) is the signal from the dc-fiber 5 with the 
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probe-tip at a specific distance from a diffuse reflecting reference material with a large, 
preferably wavelength-independent reflectance coefficient (white spectxalon) andl black 
(X) is the signal from the dc-fiber 5 with the probe-tip at that same specific distance 
from a diffuse reflecting reference material with a small, preferably wavelength- 
independent reflectance coefficient (black spectralon). Furthermore, J (X) is the signal 
trom the c-fiber 6 in contact with the sample 1 and J whU eMa& W " 




rom the c- 



fiber 6 with the probe-tip at the previously mentioned specific distance from the 
whMblack ^^atoffiTKnallyrc is a calibration constant that depends on the distance 
between the probe-tip and the reference materials. 
lQ0?Q 10 According to the invention, the processor 9 is arranged to calculate the physical 

feature using a predefined mathematicahnodel, the differential backscatter signal 
and a curve fitting mechanism. In atfeinbodiment, the diameter of the fibers 5, 6 are 
selected depending on a mean free path (mfp) of photons sent into the sample 1. It is 
noted that if the mean free path cai^tot be estimated before selecting a fiber diameter, 
15 initially two arbitrary fiber diameters may be selected. After curve fitting the measuring 
results using two different mathematical models, it will show which model amlies 
L 0°j I J Eigwes 2a and 2b depict fiber tips of the dc-fiber 5 and the c-fiber 6 in the 

situation wherein the mean free path (mfp) of photons coming out of the dc-fiber 5Kis 
much larger than a diameter d fiber of die fibers 5, 6. In anSemBoolment, the diameters of 
20 both fibers 5, 6 are of equal size^ however it should be understood that other selections 
are possible. In figure 2a, lines 21 and 22 show an example of a path traveled by a 
detectable singly scattered photon. In fig^e 2b, lines 23, 24, 25 and lines 23, 24, 26 
show two possible paths of detectable multiply scattered photons. All multiple 
scattering events occur at such large distances from the fiber tip of the fibers 5, 61that 
25 the probability of detection of multiply scattered photons is roughly equal for the dc- 
fiber 5 and the c-fiber 6. The differential backscatter signal R bs (X) will now purely be 
determined by singly scattered photons. 

f-0 0 3o J In anYembodiment, the diameteivof the fibe 



ber- 



30 In the predefined mathematical model of this embodiment, the differential backscatter 
signal R hs (X) is an exponential function of two times the mean free path. Below, an 
explanation for this model is given. 



1 00393 



7 

In the absence of absorbers, the differential backscatter signal Rbs(X) is 
proportional to the local, superficial scattering coefficient Msffl^Qsoaftyp ' As : 



Rts(X)= Cap P • -L • JdO-M^^-e^W-p-As ( 2 ) 

o 



where C app is an apparatus constant that depends amongst others on the distance 
between the probe tip and the reference materials (black and white spectralon) 5j p(^Q) 
is a function called the phase function where Q is the scattering angle, Q sca 0^& 
scattering efficiencyTpVfhe concentration of substances present in the sample 1, and As f$ 
10 the area of a scattering particle. For example, using a fused silica fiber with numerical 
aperture NA= 0.22, the differential backscatter signal R bs (X) can be approximated by 



j lit 180 

RbsW* 0^ - — - Jdcp • Jd9 ■ sin(9) . p(A 9 l 80) . ix s (A) = V p& 1 80)7*0) (3) 



0 170 



15 where q> is the azimuthal angle and 9 is the polar angle. 



dotsfofa dilute suspension of 0.2 um polystyrene spheres along with a calculation ("see 
curve 32) according to Eq. (3). In figure 3 R bs (X) is shown using arbitrary units (a.u.). 
20 Also a value of Q radar where Q radar = An- p& 1 &0yQ sca (A), is indicated in the figure. 
Ftgafe 3 shows excellent agreement between the measurement (i.e. the dots) and the 
calculation, which indicates that if mfp>d flber , the single scattering is indeed the 
dominant contributor to the differential backscatter signal R^fty as defined in Eq. (1). 
A singly scattered photon first travels from the tip of the dc-fiber 5 to a particle, 
25 and then (the same distance) hack from the particle to the tip of the dc-fiber 5 (or tip of 
c-fiber 6), see also figp-e 2a. So an average path length r(X) traveled by the measured 
single scattered photons is equal to two times the mean free path mfp(X), i.e. 



30 



z(X) = 2 • mfp#) 



(4) 



-wo immmm- c p oT fl *L20o4 /ooo 6 57 

£jQ0 H 2-3 presence of n absorbing species with specific absorption coefficients ttf^fl), 

the differential backscatter signal becomes 



n 



Rbtf) = Capp' • pa 1 80) • f4\) • exp(- <A) • X Pi • Ma**'(to) 



= Capp' • p(l 1 80) • /4(A) • exp(-2 • mfp(X)) • £ Pi • /^ ec ''(A)) (5) 



*=1 



where Capp' is an apparatus constant, p(£ 1 80) is the phase function, //y(A) is the 
scattering coefficient of the medium, X is the wavelength of the first and second 
backscattered radiation, m£p(X) is the mean free path as a function of the wavelength, n 
10 is the number of substances in the sample 1 , pi is the concentration of absorber i present 

in a detection volume of the sample 1 , and ju a spec,t Q') IS the absorption coefficient of 
absorber i as a function of the wavelength. 

It is noted that in Eq. (5) the assumption is made that absorbers are 
homogeneously distributed and do not influence each other. The Eq. (5) may be 
1 5 corrected for non-linear phenomerasuch-asaninhomogeneous distribution of 



absorbers, 

L 0 ° V V 3 According to an^embgdiment, the specific absorption coefficients of the '\^^^z^6^^~ 

absorbers, the wavelengti^dependency of the scattering coefficient /is and the phase ^ — IjJ^ 
function /?, together with Eq. (5^.are used in order to calculate the concentrations of all 
20 the absorbing substances present in the detection volume of the sample 1 . Since the 
detection volume is typically very small in the present invention, the extracted 
concentrations are highly spatially resolved. This is not possible with the known 
methods that are based on diffuse reflectance, and wherein the obtained concentrations 
are averages over large sample volumes, see e.g. Jjgp Qo 0 fetf < «/" a / < cAf<i fafffi* r 
£ QCi c/ 25 The apparatus constant Capp 9 (Eq. 3) can be determined for a specific distance 

between the tip of the dc-fiber 5 and the reference materials (black and white 
spectralon). For a suspension of monodisperse polystyrene spheres of known size and 
concentration, the scattering coefficient fi s and the phase function p(180) can be 
calculated using Mie theory^]. The apparatus constant C app ' simply follows from Eq. 
30 (3). In terms o£the volume fraction / of the suspension, the radius of the spheres a and 
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the radar efficiency coefiBcient Q ra da$)= 4tt p(X ASOyQscJX) the apparatus constant is 
determined by 



Rtfr)» Ca Pp '- J p(^180>/4(X)= Capp'- 0.05968 • Qradarfr) 

a 



(6) 
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According to another embodiment, the selected diameter dfiber is chosen so mat 
the mean free path is smaller than d/uer. In this embodiment, the differential backscatter 
signal Rbs is a function of the fiber diameter dfib er . This will be discussed in more detail 
below. 

When the mean free palh of the photons is smaller than the selected fiber 
diameter (i.e. m£p(X)<dfiber), the contribution of multiply scattered light to the 
differential backscatter signal Rb s (A) of the single dc-fiber 5 cannot completely be 
removed using Eq. (1). In this case, it appears that the average path length of the 
photons contributing to the signal Rb/A) becomes nearly independent of the optical 
properties of the sample 1. In this situation, multiple scattering events already occur at 
small distances from the tip of the dc-fiber 5. An analytical expression for the 
backscatter signal Rb s (A) is not ayailable for this situation and Monte Carlo simulations 
were used to model the beharieur ofR bs (A) as a function of the diameter of the fibers 5, 
6 and of the optical properties of the sample 1. Kgase 4 shows the results of Monte 




Carlo simulations using the MCML-code (Monte Carlo for Multi-Layered media) e# 
Wat^O-tii'[6^] for a homogeneous medium with an anisotropy value g=0.9. A flat 
circular incident beam with diameter dfiber is directed onto the sample 1, and the 
differential backscatter signal Rb S is calculated by subtracting the total reflectance in the 
c-fiber 6 (with diameter dfiber and center located at a distance dfiber from the center of the 
incident beam) from the total reflectance in the dc-fiber 5 (with diameter djtber 
overlapping the incident beam). Simulations were performed for sets of four different 
scattering coefficients (fis=15, 25, 50 and 80 mm" 1 ), four different fiber diameters 
(djiber=200, 400, 600 and 800 jam) and five different absorption coefficients (fXa=0, 0.2, 
0.4, 0.6 and 0.8 mm" 1 ). 

3 £ o o Mfr3 4 shaws R bs as a function of absorption coefiBcient |i a where the open 
circles/dashed lines correspond to ^^=200 |jm, the filled circles/dotted lines 
correspond to dfiber^OO pm, the open squares/solid lines corresp ondto j gtog600. fflgi^^ 7 ■ , # 
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and the filled squares/dash-dotted lines correspond to ^ er =800 jam. The differential 
backscatter signal Rb S for each scattering coefficient was normalized to the 
(^er=200 nm, Ha=0 mm ) case. FigBre 4 shows that in the absence of absorption, i.e. 
Pa = 0, the dififerential backscatter signal R bs depends linearly on the scattering 
5 coefficient Furthermore, the slope of the straight lines (signifying the relation 
between R bs and ^ depends only on the fiber diameter and is independent of the 
scattering coefficient The latter is more clearly demonstrated inJigjire 5, where the 
differential backscatter signal R bs is normalized to unity at zero absorption for all 
scattering coefficients |i s . The open circles correspond to 4^=200 jam and the filled 
1 0 squares correspond to 4^=800 pm. These Monte Carlo simulations therefore suggest 
that in the situation where mfp<dfib er , the diameter of the fibers 5, 6 determines the 
average path length x of the measured photons. The backscatter signal R bs for this range 
of parameters can thus be written as 



15 R bs (Z) = Cv Hs ■ exp(-T • Ha) = Cr |is • exp(-C 2 • d fiber • jUa) (7) 



where Ci and C2 are constants, x is the average path length, jia is the absorption 

coefficient, m is the scattering coefficient and d fi ber is the fiber diameter of the fibers 5, 
6. 

V00 ^ff* CXaCt analytica ^ ex P ress i° n for R b$) is not available due to the large 

contribution of multiple scattering events to the signal. Measurements were done for 
determining a total integrated backscatter signal R iot for a range of X between 400-900 
nm, using the formula 



900 mn 



25 Rtot(Ms)= f dXRtAMs) (8) 



400nm 



£ Q 0J7^23 6 shows ^ mte grated total backscatter signal R to t(tt) is proportional to 

fi s (X) in the relevant range of & between 10-100 mm" 1 . Therefore, in the absence of 
absorbers, it follows that 
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(9) 



which is in agreement with the Monte Carlo simulations. 

In the presence of n absorbing substances in a suspension, with specific 
absorption coefficients fia^fX), the differential backscatter signal becomes 



n 



(10) 



10 



where r is the average path length of the detected backscattered photons and p f is the 
concentration of the substance i. 
POO f?-^} Non-linear phenomena such as an inhomogeneous distribution of absorbers are 



15 



not incorporated in Eq. (10), but can be added by the skilled person, see e.g. VtM 

Figure 7 shows the measured and calculated average path length r as a function 
of the average scattering coefficient <Hs(k)> (with 500 ra 




mil 



m < X < 700 nm) for d fiber = 0.4 
_1 at A=600 nm. In Seare 7 



20 



25 



30 



and for absorption coefficient jUafi) = 2.0 m 
measurement are indicated by dots, the r = 2-mfp curve is indicated by a line 71 and 
Monte Carlo simulations are depicted by dashed lines. Identical results were obtained 
for suspensions with an absorption coefficient of = 1 .0 mm ' 1 at 600 nm. The average 
path length r was determined using suspensions of polystyrene spheres with different 
sizes and concentrations to vary the scattering coefficient The anisotropy g of 
these suspensions was in the range of 0.8-0.9. Evans Blue dye was added as an 
absorber, and the average path length r was calculated from Eqs. (9) and (10) and 
knowledge of the concentrations an^pecific absorption coefficient of Evans Blue, as 
will be known to the skilled person. 

Looking at the measured average path lengths of figase 
large scattering coefficients (//, = 10-100 mm" 1 , the range relevant for tissue) the 
average path length r is independent of the scattering coefficient fi s to within 10% and 
approximately equal to half the fiber diameter (r * 0.24 mm while d fiber = 0.40 mm), in 
agreement with the Monte Carlo simulations (the dashed lines correspond to Monte 
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Carlo calculations for d fibe r = 0.2, 0.4, 0.6 and 0.8 mm). For small scattering coefficients 
(e.g. /4 < 5 mm" 1 ) the average path length r is well described by r= 2-mfp, according 
to Eq.(4), see line 71 . Kguie 7 also clearly demonstrates that the transition from the 
'single scattering regime' to the 'constant path length regime' occurs for mean free 
5 paths of the order of the fiber diameter. It is therefore expected that single scattering 
will prevail over a larger range of scattering coefficients for fiber diameters smaller 
than. 400 pm. 

In the following, the effect of absorption on the average path length r will been 
examined in more detail. Various concentrations of Evans Blue dye were added to a 
10 suspension of polystyrene spheres with scattering coefficients /4 of 35 mm" 1 . The 

concentrations of Evans Blue (EB) dye were varied such that the absorption coefficient 
fia at 600 nm was in the range of 0 to 2 mm" 1 . Typical results of the differential 
backscatter signal R bs for three different absorption coefficients are shown in figare ^ 
8. Note that the signal I of the dc-fiber 5 is plotted on a different vertical scale than the 
1 5 si^l J of the c-fiber 6 and the differential backscatter signal R bs . 

ipectra with Evans Blue R EB present in the suspension were divided by the 
spectrum with no Evans Blue present R° and the negative natural logarithm of the ratio 
R EB /R° was determined: 




20 A = - In^/R 0 ) = x • p • fj^ cMJt (1 1) 

where p is the concentration of Evans Blue, and jUa spe<iEB is the specific absorption 
coefficient of Evans Blue. 

\ 0p5 7\ 25 Ijigare 9 shows a typical spectrum of an absorption curve 92, along with the 

specific absorption coefficient /fc*"** 8 of Evans Blue dye, see curve 94. 

For all concentrations, an area A* under the absorption curve 92 was determined 
in the wavelength range X between 500 and 650 nm. From Eq. (1 1) it follows that if the 
average path length x is independent of the absorption coefficient area ,4* should 
30 depend linearly on the concentration p of the species in the suspension. 
£ (y>T9 3 ^gtaSe 1 0 shows a measured A * as a function of the absorption coefficient ^ at 

600 nm. A curve fitting is done, resulting in a line 104 for ^(X )=35 mm' 1 . 10 
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shows that the average path length x is indeed independent of the absorption coefficient 
Ha in the range 0-2 — "* 



iiiii 



£ QOvjP J From the previous results of figure 4 to 1 0 5 it shows that for mfp < ^ e r> the 
differential backscatter signal R bs is described by Eq. (7) with C 2 * 0.6. Figujf£ 1 1 
5 shows typical spectra measured in a suspension of 1 .0 jim polystyrene spheres with and 
without Evans Blue dye 2 and 0 mm" 1 at 600 nm, respectively). From Eq. (7) and 
.figure 7, the relation between the differential backscatter signals with and without 
absorber is given by 



10 Rbs(Kf^} = R bs (^0yexp(-0.24^ (12) 

qQ 6 ( J The calculated spectrum according to Eq. (12) is plotted as a dashed line 1 10 in 
fij*Bre 1 1 and shows excellent agreement with the measured R bs (A,/Ja), see line 1 1 1 . In 
figHse 1 1, line 1 12 depicts R bs & 0). 
15 In short, the average path length t of photons measured when subtracting the 

^ signals of the c-fiber 6 from the dc-fiber 5 using Eq. (l)Vis independent of the optical 

properties of the sample 1 and approximately equal to naif the diameter of the fibers 5, 
6 used, as long as the fiber diameter is larger than the mfp. 

20 In a specific embodiment, the device according to the invention is arranged to 

)^Oo(q ^3 determine concentrations of oxygenated blood in tissue. Since the scattering coefficient 

of tissue juj issue is in the range of 10-100 mm" 1 , the fiber diameter should be smaller 
than a certain maximum diameter dmax where dmax is between 10 and 100 |om 9 ^for 
example^smaller than 50 nm, in order to measure predominantly single scattering in 

25 tissue. In this case, Eq. (5) holds. For fibers 5, 6 with much larger diameters (e.g. 200 
or 400 nm), the differential backscatter signal R bs (X) is described by Eq. (10) with r« 

j^O 0 G ^3 * S present * y k* 10 ™ 11 th at th e wavelength dependence of the scattering coefficient 

in tissue jUg ilsslte can be adequately described by an empirical power-law function, see______^ 

71* Hh p(e M${ oAcle ctkJ[ hctfo, >j A Ayp^L i] o» 
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with a and b constants that depend on the size, concentration and relative refractive 
index of the scatterers (i.e. substances) present in the detection volume. 
f {\ft(o $3 [iSfdominant absorbers in tissue in the visible wavelength range are oxygenated and 
5 deoxygenated blood. Thus in tissue Eq. (10) becomes 

Rt&) = Capp • aT b • exp(- • Pblood - (S 02 • ^ ox + (1-S 02 ) • a*"**" )) 



= CVp ■ ■ exp(- 0.6-djuvr • Pbiood • (S02 - fiT** + (I-S02) • ^ ec4eox )) (14) 



10 



where pwood is the concentration of blood, S02 is the blood oxygenation (percentage 
oxygen saturation) in a c^tain detection volume, 6app is a constant that depends on the 
calibration constant c, C£a$p is C app -a, X is the wavelength, b is the slope of the 
scattering coefficient defined in Eq. 13, fjj pec,0x is the specific absorption coefficients of 
15 fully oxygenated b\oody^ spec,deox is the specific absorption coefficients of fully 
deoxygenated blood. 

£ 0 0 (2(9 3 Non-linear phenomena such as an inhomogeneous distribution of absorbers are 

not incorporated in Eq. (14), but can be added by the skilled person, see e.g. [8]. 
\ O & G 7 3 Since the specific absorption coefficients of fully oxygenated (jUa^ ec,0X ) and fully 

20 deoxygenated (jUa 5pe *) blood are well known, see Sgdre 12, Eq. (14) can be fitted to 
the measured data to yield the slope b of the scattering coefficient jUs tissue , the 
concentration pbiood and the oxygen saturation S02 of the blood present in the detection 
volume. When a correction is made for the inhomogeneous distribution of blood in the 
vessels, a vessel diameter D may be determined as well. Since the average detection 
25 depth is small (e.g. 0. 1 mm), the blood present in the detection volume when measuring 
non-invasively is located in capillaries. 
1 0 Cio^ In fijft 13, in vivo measurements of backscattering in ■ a human trachea together 

with a fit using Eq. (14) are shown. The measurements are realized using a fiber 
diameter of 400 pm. The dots depict the measurements and curve 130 is a fitting curve. 
30 In fi^e 1 3, b - -0.94 and the oxygenation S02 = 95%. 

The present invention can be used for tumor detection. Tumor growth may, due 
to its excessive oxygen consumption, be accompanied by a low capillary oxygen 
saturation, which can only be assessed using a very localized measurement. Since (pre- 
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)cancerous tissue is generally more heterogeneous than normal tissue, the standard 
deviation of multiple measurements is likely to be larger for (precancerous tissue than 
for normal tissue. Standard deviations in the measurements can be calculated for the 
oxygen saturation, the blood concentration, the blood vessel diameter, and the slope b 
5 of the scattering coefficient It is noted that the invention is by no means 
restricted to determine a concentration of a substance as the physical feature. All 
features, mentioned in the previous phrase can be regarded as physical features. 



An example of a measurement of a lung tumor is shown in^Sgwe 14. The shape 
10 of the dip in the wavelength range of 500-600 nm in this figure demonstrates the 
depletion of oxygen from the capillaries of this tumor due to its excessive oxygen 
consumption. 

0 7 0 When a needle-probe is used, the local oxygenation and scattering coefficient 

fis tissue can be measured invasively. This could be helpful in determining tumor-margins 
15 intra-operatively in real-time, for instance during resection of a breast-tumor. 
£ Ofl 7 ^ J According to an embodiment, the device comprises multiple probes and a 
multichannel spectrometer for multiple simultaneous measurements on different 
locations of the sample 1. Using this device, multiple measurements can be made 
— simultaneously on different locations of for example a suspicious lesion. 

20 — ; In yet another embodiment, the device comprises at least two pairs of fibers, 

having different fiber diameters. For example, when a pair of fibers with 100 jxm, a pair 
of fibers with a diameter of 200 nm and a pair of fibers with a diameter of 400 pm are 
used, information from different depths in the sample 1 can be obtained as the average 
_ path length increases with increasing fiber diameter. 

25 J The method and apparatus according to the invention can also be used to analyze 
local drug concentrations. From Eq. (10) it follows that if the specific absorption 
coefficient of a certain drug is known, the local concentration p of that substance can 
be determined using the invention. 
£q A *7Sji Another possibility of the present invention is to monitor glucose concentrations. 
30 The scattering coefficient pf™" depends among others on the relative refractive index 
of the scatterers with respect to the surrounding medium (in tissue: cytoplasm). The 
refractive index of the surrounding cytoplasm is likely to depend on the concentration 
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of glucose. A change in the glucose concentration will therefore likely affect the slope 
b of the scattering coefficient p/""* , see Eq. (13). 
] 0 y(p J Wfcfet specific embodiments of the invention have been described above, it will 
be appreciated that the invention may be practiced otherwise than as described. For 
5 example, a concentration of a substance in polluted water may be calculated. The 
description is not intended to limit the scope of the invention. 



10 



9051 



17 



1. A. AmeluJc, M.P.L. Bard, J A. Burgers and H.J.C.M. Sterenborg, "Single scattering 
spectroscopy for the endoscopic analysis of particle size in superficial layers of turbid 

5 media", Apphed^Optics 42, page 4095-41 0 1 (2003). 

2. R.M.P. Doornbok R Lang, M.C. Aalders, F.W.Cross and HJ.C.M. Sterenborg, "The 
determination of in viV) human tissue optical properties and absolute chromophore 
concentrations using sopfcrtially resolved steady-state diffuse reflectance spectroscopy", 
Phys. Med. Biol. 44 (1999>oage 967-981. 

10 3. A. Angstrom, "On the atmospheric transmission of sun radiation and on dust in the 
air", Geograf. Ann. Deut.ll, 15^- 166 (1929) 

4. H.C. van de Hulst, "Light Scattering by Small Particles", Wiley, New York, 1957. 

5. R. GraafF, J.G. Aarnoudse, J.R. ZijV P.M.A. Sloot, F.F.M. de Mul, J. Greve, and 
M.H. Koelink, "Reduced light-scatteringjproperties for mixtures of spherical particles: 

15 a simple approximation derived from Mie calculations", Applied Optics. 31, 1370-1376 
(1992). 

6. L-H Wang, S.L. Jacques, L-Q Zheng, "MCML - Monte Carlo modeling of photon 
transport in multi-layered tissues", Computer Methods and Programs in Biomedicine 
47, 131-146 (1995). 

20 7. L-H Wang, S.L. Jacques, L-Q Zheng, "CONV - CoiWolution for responses to a 

finite diameter photon beam incident on multi-layered tis^jxes", Computer Methods and 
Programs in Biomedicine 54, 141-150 (1997). 

8. R.L.P. van Veen, W. Verkruysse, and HJ.C.M. SterenborgVDiffuse reflectance 
spectroscopy from 500 to 1 060 nm using correction for inhomogbneously distributed 
25 absorbers", Opt Lett. 27, 246-248 (2002). 



£LQ-2Mm29QSl___ E €T/NL20Q4/000657 - 

1 . A method of determining of a physical feature of a medium, comprising: 

- producing radiation with a light source (2); 

5 - placing a probe on a sample (1) of said medium, said probe comprising a first optical 

\ 

fiber (5) having a first diameter, and at least a second optical fiber (6) having a second 
diameter; 

- sending light coming from said light source, through said first optical fiber; 

- collecting first backscattered radiation through said first optical fiber and second 
1 0 backscattered radiation through said second optical fiber; 

- producing a first signal^ (I) based on said first backscattered radiation, and a second 
signal (J) based on said second backscattered radiation; 

- determine a measured differential backscatter signal as a function of wavelength using 
said first and second signals y, J), 

15 characterized by 

- calculating said physical featur^by curve fitting said measured differential 
backscatter signal to a backscatter>function, in which said backscatter function is a 
function of an average path-length (x) travelled by detected scattered photons, said 
average path-length (x) being independent from an absorption coefficient (jig) of said 

20 medium, and from a scattering coefficient (jig) of said medium. 

2. Method according to claim 1 , wherd^said average path-length (x) is also 
independent from a wavelength (X) of said firstand second backscattered radiation 




25 3. Method according to claim 1, wherein said path-length (x) is proportional to said 
first fiber diameter. 



4. Method according to claim 1, wherein said backscatter function is given by: 
30 Rbs^l-jHs-exp^x/fia) 

With X = C2-dfiber 
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where CI and C2 are constants, fi a - said absorption coefficient of said medium, ^ = 
said scattering coefficient of said medium, and dfiber = said first fiber diameter. 



5 



5. Method according to claim 4, wherein C2 is approximately 0.6. 



6. Method according to any of claims 1-5, wherein said physical feature is a 
concentration of at least one substance in said medium. 



7. A device for determining a physical feature of a medium, comprising: 
10 - a light source (2) for producing radiation; 

- a probe with a. ■J^U and a second optical fiber (5, 6), said firs, optical fiber (5) 
having a first diameter and being arranged to deliver said radiation on a sample (1) of 
said medium and to collect first backscattered radiation from said sample (1), said 
second optical fiber (6) haying a second diameter and being arranged to collect second 

15 backscattered radiation, wherein said second optical fiber (6) is positioned alongside 
said first optical fiber (5); 

- a spectrometer (7) for producing a first signal (I) based on said first backscattered 
radiation, and for producing a second signal (J) based on said second backscattered 
radiation; 

20 - a processor (9) arranged to determine a measured differential backscatter signal as a 
function of wavelength (A) using saifcl first and second signals (I, J), 
characterized in that said processor i^arranged to calculate said physical feature by 



curve fitting said measured differentiaTbackscatter signal to a backscatter function 
(Rbs), in which said backscatter functionis a function of an average path-length (x) 
25 travelled by detected scattered photons, said average path-length (t) being independent 
from an absorption coefficient (jia) of said medium, and from a scattering coefficient 
(jis) of said medium. 



8. Computer program product to be loaded by a^computer, said computer program 
30 product, after being loaded, providing said computer with the capacity to: 

- receive a first signal (I) indicative of a collected radiation received from a first fiber 
(5) and a second signal (J) indicative of a collected radiation received from a second 
fiber (6); 
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- determine a measured differential backscatter signal (R^) as a function of wavelength 
(A) of said collected radiation using said first and second signals (I, J); 
characterized\bv the capacity to 

- calculate a physical feature by curve fitting said measured differential backscatter 
5 signal to a backscatter function, in which said backscatter function is a function of an 

average path-length (x) travelled by detected scattered photons, said average path- 
length (x) being independent from an absorption coefficient (p-a) of said medium, and 
from a scattering coefficient ([is) of said medium. 



10 9. Data carrier provided with a computer program product according to claim 8. 



10. A method of determining ^physical feature of a medium, comprising: 

- producing radiation with a light source (2); 

- placing a probe on a sample (1) of said medium, said probe comprising a first optical 
15 fiber (5) having a first diameter, and at least a second optical fiber (6) having a second 

diameter; 

- sending light coming from said light source, through said first optical fiber; 

- collecting first backscattered radiation through said first optical fiber and second 
backscattered radiation through said second\optical fiber; 

20 - producing a first signal (I) based on said first backscattered radiation, and a second 
signal (J) based on said second backscattered radiation; 

- determine a measured differential backscatter signal as a function of wavelength using 
said first and second signals (I, J), 
characterized bv 

25 - calculating said physical feature by curve fitting said measured differential 

backscatter signal to a backscatter function, in which said backscatter function is a 
function of a mean free path of photons. 



1 1 . Method according to claim 1 0, wherein said backscatter function (Rbs) is defined 
30 by: \ 

Rtefl.) = Capp' • 180) • Hs(X) • exp(-2 • mfpQ)) ■ £ Pi • ^ spec,i (X)) 
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where W is aa apparatos cnston., *J 80) is a phase Morion, is a scattering 
coefficient of said medium, X is a wavelength of said first and second backscattered 
radiation, mfp(A.}is said mean free path as a function of the wavelength, n is a number 
of substances in said medium, pi is concentration of absorber i present in a detection 
5 volume of said sample (1), and Ha Spec,1 (X,) is an absorption coefficient of substance i as a 
function of the wavelength. 

12. Method according^ any of claims 10-11, wherein said physical feature is a 
concentration of at least one substance in said medium. 

10 

13. A device for determining a physical feature of a medium, comprising: 

- a light source (2) for producing radiation; 

- a probe with at least a first andV second optical fiber (5, 6), said first optical fiber (5) 
having a first diameter and being arranged to deliver said radiation on a sample (1) of 

15 said medium and to collect first backscattered radiation from said sample (1), said 

second optical fiber (6) having a second diameter and being arranged to collect second 
backscattered radiation, wherein said second optical fiber (6) is positioned alongside 
said first optical fiber (5); 

- a spectrometer (7) for producing a first signal (I) based on said first backscattered 
20 radiation, and for producing a second signal v (J) based on said second backscattered 

radiation; 

- a processor (9) arranged to determine a measured differential backscatter signal as a 
function of wavelength (X) using said first and second signals (I, J), 

characterized in that said processor is arranged to calculate said physical feature by 

\ 

25 curve fitting said measured differential backscatter signal to a backscatter function 
(Rb S ), wherein said backscatter function is a function ofa mean free path of photons. 

14. Computer program product to be loaded by a computer, said computer program 
product, after being loaded, providing said computer with the capacity to: 

30 - receive a first signal (I) indicative for a collected radiation received from a first fiber 
(5) and a second signal (J) indicative for a collected radiation received from a second 
fiber (6); \ 
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- detennine a measured differential backscatter signal (Rbs) as a function of wavelength 
(k) of said collected radiation using said first and second signals (I, J); 
characterized by the capacity to 

- calculate a physical feature by curve fitting said measured differential backscatter 
5 signal to a backscatter function, wherein said backscatter function is a function of a 

mean free path of photons. 



1 5 . Data carrier provided with a computer program product according to claim 1 4. 

10 16. Method according to any of the claim 1-6, 10-12 wherein said method comprises: 

- simultaneously measuring backscatter radiation on different locations of said sample 

(i); 

- determining a physical feature for said different locations; 

- calculating a standard deviation of said physical feature. 

15 
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